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Abstract 

We are concerned with the problem of finding among all polynomials of degree at most 
n and normalized to be 1 at c the one with minimal uniform norm on £. Here, £ is 
a given ellipse with both foci on the real axis and c is a given real point not contained 
in £. Problems of this type arise in certain iterative matrix computations, and, in this 
context, it is generally believed and widely referenced that suitably normalized Chebyshev 
polynomials are optimal for such constrained approximation problems. In this note, we 
show that this is not true in general. Moreover, we derive sufficient conditions which 
guarantee that Chebyshev polynomials are optimal. Also, some numerical examples are 
presented. 
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1. Introduction and Statement of the Main Results 


Let U n be the set of all complex polynomials of degree at most n. For r > 1, we denote by 

£ r := {z £ C | \z — 1| + \z + 1 1 < r -} — } 

r 

the ellipse with foci at ±1 and semi-axes 

a r + ) j ^( r ) 

l r 2 r 

In this note, we study the constrained Chebyshev approximation problem 


% 


min 

pen n :p(c)=i 


max \p(z)\ 


(i) 


where n £ IN, r > 1, and c<ElR\£ r . Standard results from approximation theory (see e.g. 
[9]) show that there always exists a unique optimal polynomial, denoted by p n (z;r,c ) in 
the sequel, for (1) and, moreover, that p n is a real polynomial. In 1963, Clayton [3] proved 
that p n (z;r, c) is just the polynomial 


t n (z-,c) 


Tn(z) 

Tn(c) 


( 2 ) 


where 

T.<*) = ;(<’" + 4> • * = 5(0 + ^) ■ (3) 

2 v n 2 v 

denotes the nth Chebyshev polynomial. The approximation problem (1) arises in certain 
iterative matrix computations (see e.g. [2,5]). In this context, Clayton’s result is widely 
referenced in the literature (e.g. [2,5,8,12,13]) and is even used to derive new results 

on constrained approximation problems [1]. Surprisingly, nobody seems to have checked 
Clayton’s proof. 

In this note, we show that the normalized Chebyshev polynomials (2) are not always 
optimal for (1), and hence Clayton’s result is not true in general. More precisely, we have 
the following 


Theorem 1. 

a) Let r > 1 and c > a r or c < —a r . Then , for n = 1,2, 3, 4, t n (z]c) is the unique 
optimal polynomial for (1). 

b) For any integer n > 5 there exists a real number r* = r*(n) > 1 such that t n (z;c) is 
not optimal for (1 ) for all r > r* and all c € 1R with a r < \c\ < a r + l/&r* 

However, t n = p n in most cases, and t n ceases to be optimal only for normalization points 
c which are very close to the ellipse. We will show that the following conditions on c are 
sufficient to guarantee the optimality of t n . 
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Theorem 2. Let n > 5 be an integer, r > 1, and c <E JR. Then, t n (z;c) is the unique 
optimal polynomial for (1), if 

(a) |c|>i(r^ + r-^) 
or 

( b) \c\ > (2«r - 1 + \/2«r - a? + 1) . 

Remark 1. In general, the conditions (a) and (b) do not imply each other. In particular, 
(a) (resp. (b)) is less stringent for small r (resp. large r). Also, note that (b) is satisfied 
if |c| > (1 + y/2/2 )a r . 

The paper is organized as follows. In Section 2, we state a necessary and sufficient 
criterion for t n to be optimal for (1). Also some auxiliary results are collected which will 
be used in Section 3 and 4 to prove Theorem 1 and 2, respectively. Finally, in Section 5, 
we present some numerical examples. 

2. Preliminaries 

In the sequel, let always be r > 1 and n € IN. Since p n { z > r i ~ c ) = Pn {~ z j T \ c ) h * s sufficient 
to consider positive c only; so for the rest of the paper, wo assume that c > a, . 

First, we determine the extremal points zi of t n defined by 

\t n (zp,c)\ = max |<„0;c)| , z i € £ r ■ 

With (3), one easily verifies that there are 2n such points given by 

Z) := a T cos <pi + ih r sin <pi , := hr/n , l = l,...,2n 

Moreover, note that t n (zf\c ) = ( — l) i T„(a r )/T n (c). Using' Rivliu and Shapiros character- 
ization [10] of the optimal solution of general linear Cliebyshev approximation problems, 
we deduce that t n = p n iff there exist nonnegative real numbers ai, l = 1, . . . , 2n (not all 
zero) such that 

2n 

V <T/(-l ) l q(z[) = 0 for all </ e IT n with q(c) = 0 . (4) 

/=i 

By solving this linear system explicitly, one arrives at the following 
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Lemma 1. The polynomial (2) t n is optimal for (1 ) iff > 0 for 2 n, where 




T fc (a r ) 




( 5 ) 


Proof. The result is a special case of Theorem 3 in [4] where we investigated the ap- 
proximation problem (1) in the more general setting of complex c. On the other hand, 
by using the polynomials q(z ) = Tk(z) — T*(c), fc = l v , M n, as a basis in (4), it is also 
straightforward to verify directly that the <t/ given by (5) satisfy (4) and that these are up 
to a constant factor the only solutions of (4). ■ 

Remark 2. Clearly <7 2n > 0 and, moreover, <j\ = cr 2 Hence, t n is optimal iff > 0 
for / = 1, . . . , n. 

The following result due to Rogosinski and Szego [11] will be used in the next section 
to establish a sufficient condition for the positivity of the <7/. 

Lemma 2. Let Ao, Ai,. . . , A„ be reai numbers which satisfy A n > 0, A n _] — 2A n > 0, and 
Afc_i — 2A* + Ajt-h i ^ 0 for k = 1 , 2, , . . , n — 1 . Then: 


*(*>) ■■= t + E A* cos(kif) > 0 for all <p E 1R 


( 6 ) 


k=l 


We close this section with the following technical lemma. The proof is straightforward 
and omitted here. 


Lemma 3. 

a) Let k E IN- Then: 

f «o t V-W = 1° i{ k = 1 

k \ k/2 if k> 2 

j=i v ' — 

b) Let 2 < / < n be an even integer and < pi = ln/n. Then: 


and 


n — 1 

^ ^ cos(fci^i) = 0 

k=0 

TJ — 1 

y kcos(k<pi) = — n/2 

Jt=i 


( 7 ) 

( 8 ) 
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3. Proof of Theorem 1 


Let r > 1 be fixed and set a := a r . Then, for each /, (5) defines a polynomial <r ( (c) = cr, 
in c of degree n. Therefore, 

71 ( j) , \ 

<ri(c) = <r t (a) + (c- a)(a\(a) + "4 — ( c - a ) 7 ') • ( 9 ) 

i — 2 ,/ ' 

First, we prove part b) of Theorem 1. Let n > 5 and 2 < / < n be an even integer. With 
(5) and (7), it follows that 

71—1 

a, (a) = (-1)' (- (1 + (-!)') + Y cos ( fc< ^)) = 0 • ( 10 ) 


Furthermore, we derive from (5) 


a',{a) 


I + V 

2 T n(a) rrj 


Tt(«) 


cos(ktpt) . 


( 11 ) 


Let = cos((2j - 1 )tt/(2As)), j = 1, . . . , fc, denote the zeros of T*. Then, 


*( a ) _ v' i _ V" 1 Vb w r V - Vff (fc) ) 2m 

\7T\ ~ 2_> c(k) ~ 2L> a m+l a 2m+J Z-Z^J ' 

k{a) +r[ a-Cj m=0 a J=1 m= o J=1 

0 if k = 1 

ib/(2a 3 ) + 0(l/a 5 ) if A: >2 

Here, we used the fact that T(./Tfc is an odd function and part a) of Lemma 3. With (8), 
(11), and (12), it follows that 


= k/a + 


( 12 ) 




(13) 


Combining (9), (10), and (13) yields 

1 /7T , 1 


" (j) 


- (c) =(c -„)(-! co4)i + 0(i, + E^ 

j — 2 

and, finally, since, in view of (5) and T^(a)/Tr(a) = 0(l/a.i), f or j > 2 we have <r^(a) 
<9(l/u 2 ), 

*,w = ^ (-r os O + 0( ^ ) + ° Wc _ a,) ) ■ 
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Thus, 07 (c) < 0 and, therefore, ( 2 ) is not the optimal polynomial for (l), if c — a < 1 /a 2 , 
a is sufficiently large, and cos(hr/n) > 0 , i.e. I < n/2. Note that even / with 2 < / < nj 2 
exist, since n > 5. This concludes the proof of part b) of Theorem 1 . 

We now turn to the proof of part a) of Theorem 1 . Let r > 1 and c > a = a r be fixed. 
Moreover, set A * := T*(c ) and a* := 2*(a). Then, in view of Lemma 1 and Remark 2 , one 
needs to check the positivity of 


^ = M)' (j f 1 +(- 1 >'^) + E £ cos O )• , = 1 -- n ■ 

n *=1 * 


(14) 


for the four cases n = 1 , 2 , 3, 4. For n = 1,2 this is clearly true, since 

d !) = h— - !) > 0 » d 2) = - 1 ) > o , 

2 aj i fl2 


and 


,«> = _ 2^i + 1) = (e - ° K . a r g + !) > o 

2 a 2 a] 


a( 2 a 2 - 1 ) 

Next, consider n = 3. It is easily verified that ,4 3 /a 3 > ^4j/a 1? and hence 


.(3) 


= |(— -l )>0 

2 03 Q \ 2 a 2 


By using that T 2 (c)T 2 (a) + ca is a monotonously increasing function in c for c > a > 1 , we 
deduce 


( 3 ) 1/^3 yl 2 A\ -x (c-a)/2T 2 (c)T 2 (a) + 2 ca + l 

^ 2 “ 2 l a 3 “ a 2 " a, j “ 2 a V (4a 2 - 3)(2a 2 - 1 ) 

2 a(c — a) 

> > 0 . 

“ (4a 2 - 3)(2a 2 - 1) 


Similarly, one obtains 


( 3 ) 1 .4 3 A 2 Aj 1 (c - a) / 4(c 2 + 2 ca + a 2 ) - 3 

<7s 2 a 3 a 2 ai 2 a \ 2(4a 2 — 3) 

^ (c — a)(16a 4 -18a 2 + 9) ^ 

- 2a(4a 2 - 3)(2a 2 - 1) > ° ' 


2 ca + 1 \ 
2 a 2 - 1 ) 


Finally, we turn to the case n = 4. Analogously to the case n = 3, / = 1 

2 CL 4 2 {Z 3 CL\ 
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For / = 2, we have 


(4") 1,A» (A -2 - a 7 ){A 2 a 2 - a\ + 1) ^ 

o \ — -( 2 h 1) = 777-2 77 > u 

2 a,j a2 ^2(2a 2 — 1) 


The positivity of <73^ follows from 
A*) 


— a 2 ) V a* 03 a \ ' 


2 (c 2 — a 2 ) 4(c 2 — a 2 ) V 04 

2(c 2 + a 2 — 1) \/2c 


8a 4 — 8a 2 + 1 a(4a 2 — 3) 

8(2 - y/2)a 4 + 4(2 y/2 - 5)a 2 + 6-^2 > () 


(15) 

(1C) 


(8a 4 -8a 2 -f l)(4a 2 -3) 

Here we have used that (15) is a monotonously increasing function in c for c > 1 and that 
the numerator in (16) has no real zero. Similarly, by a routine, but lenghty, computation, 
one verifies that 

a 2 a 3 a 4 ^(4) _ 02^3^4 _|_ 7^2 _ Aj_ ^ T\ 

2(c — a) °" 4 _ 2(c — a) \2 a 4 a 3 a 2 a, 2/ 

= (2c 2 — l)^(c — a)a 3 + a 2 ^a 2 + ^(c(4a 2 — 1) — «s)( a 2 ~ l) a ~ a 2 j( a 2 — 1) 

x* a 2 (4a 4 — 6a 2 H— 3 ) — (- 2a 2 (a 2 — 1) >0 

This concludes the proof of part a) of Theorem 1. 

4. Proof of Theorem 2 

Let r > 1 and c > a := a r be fixed. Note that a and c have the representations 


„ = !(,. + 1) , 0 = 1(7? + !) , /?>■• 


With (3) and (17), one obtains 

T k (c) R k + l/R k 


= f(<Pk) , 


where we set 


Tk(a ) r k + l/r k 

cosh((log R)mp/ir) _ for 

^ cosh((log r);?<^/7r) ’ < ^ >k n 


(17) 


(18) 


Since / is continuous, bounded, and even, it can be expanded into the Fourier series 

1 00 

f{<p) = -«0 COs(jV) , -7T < If < 7T . 


J=1 
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By rewriting the expression (5) for < 7 / in terms of (18) and, subsequently, using the discrete 
orthogonality relations of cos {l<Pk), k,l = 0, . . . ,n, (see e.g. [7], p.472), we get 

<n = (-i)'(km + (— 1)7W) + E SM <**(>?*)) 

k=l 

j f(-l)'(a/ + X^m=i( Q 2 mn-; + a 2 mn+()) for / = 1 , . . . , n - 1 
(«(-l)'(« n +Em=l°2(m+l)n) foi l = Tl . 

It follows that all <n > 0 and, in view of Lemma 1 , that the normalized Chebyshev 
polynomials ( 2 ) are optimal for ( 1 ), if the Fourier coeffients of / satisfy 

a } = > J — 1)2,... . (19) 

It is well known (see e.g. [ 6 ], Theorem 35) that (19) holds true if / is a convex function. 
Hence, in order to prove that the condition (a) in Theorem 2 guarantees the optimality 
of the polynomial ( 2 ) for ( 1 ), it only remains to show that (a) implies the convexity of 
/. Since / is even, we only need to consider <p > 0 . Moreover, set x (\ogr)ntp/i: and 
7 := log R / log r > 1 . Then, using standard calculus, we obtain 

j- \ ( — ) f"(<p) = 7 2 — 1 — 27 tanh(x) tanh( 7 ar) + 2 tanh 2 (a i ) 

cosh( 7 a;) Vnlogr/ 

> 7 2 — 1 — 27 tanh(x) + 2 tanh 2 (x) 

> 7 2 - 1 + 2 min y{xj - 7 ) ( 20 ) 

0< V <1 

__ f (1 - 7 ) 2 if 7 > 2 
\ 7 2 /2 — 1 if 7 < 2 . 

Therefore, ( 20 ) is nonnegative, and thus / convex, if 7 > \/2. This last condition is easily 
seen to be equivalent to the condition (a) in Theorem 2 . 

Remark 3. The main idea of the proof, namely to verify the positivity of the at via the 
convexity of /, is due to Clayton [3]. However, in [3], it is claimed that / is convex in all 
cases R > r > 1 . Unfortunately, this is not true in general. 

Now, assume that the condition (b) of Theorem 2 is fulfilled. Again, we will use the 
notations A* = 7fc(c) and a = Tk(a). Note that, by the three-term recurrence formula of 
the Chebyshev polynomials, 

Ak+i = 2cAk — Ak-i , k = 1 , 2 ,... . ( 21 ) 
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Next, set 


A 0 = — , A n = ]-, and , for k = 1, 2, . . . , n - 1 , A fc = 

a n 2 &n-k 


(22) 


and let s(ip) be the trigonometric polynomial defined by (6). With (5) and (6), one readily 
verifies that <t; = s(ln/n), and, in view of Lemma 1 and 2, we conclude that the polynomial 
(2) is indeed optimal for (1) if the numbers (22) satisfy 

A„ > 0 , — 2A n > 0 , and , for k = 1, . . . , n — 1, A^._, — 2A* + Ajt +] > 0 . (23) 

The first condition in (23) is trivially true, and the second one follows from A\ > a*. Using 
(22), the remaining inequalities in (23) can be rewritten in the form 


and 


Aj+1 _ 2 ^i + 


A simple calculation shows, that (24) is equivalent to 

a 2 + \/a 2 a 2 + 1 


— — 2— 

+ 1 >0 

(24) 

«2 a l 

2 


A ‘- x > 0 

, for j = 2, . . . , n - 1 . 

(25) 


c > c* := 


2 a 


(= 


2a 2 — 1 + y/2 a* - a 2 + 1 
2a r 


) 


(26) 


which is just the condition (b). For the proof of Theorem 2, it only remains to show that 
(26) also implies (25). Let j > 2. First, by using (21), we deduce that 

A>+i _ 2 ^ 2 . -(- A} ~ 1 

(Ij & j 2 


a j + 1 


> 


( , C 1 . 1,1 1 A A,-2 / 1 1 \ 

M 2( ^ “ 7? + 7 " + 

- (^4c 2 a,rtj_i — 4caj+iaj_i + aj(aj+i — a^-i)) 

j+l a j a j-i ' 


(27) 


2cajq_i< 

Next, set 

<?j(c) 4c 4 cOj_|_|Uj — i F — l) 

and note that Q } attains its minimum at aj+i/(2aj) < c* . Hence, in view of (27), (25) 
holds true, if Q } (c*) > 0 is fulfilled. This is indeed the case, and we will show by induction 
that 

Q 3 (c) > Q 2 (c*) > 0 , j = 2,3,... . (28) 
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For j = 2 , this follows with 

Q 2 (c*) = 4(c*) 2 a 2 a — 4 c*a 3 a + a 2 (a 3 — a) 

= a -1 ^a 2 (2a 4 — 3a 2 + 2) — ( a 2 — l)\/a^a 2 + 1^ > 0 , 

since y/2a 2 > \/a^a 2 + 1 and 2 a 4 — 3 a 2 + 2 > \/2(a 2 — 1 ) for a > 1 . Finally, if (28) holds 
true for j , a routine, but lengthy, calculation shows that 

<?>+i(c*) - Qi(c‘) = (<H - l)(-4(c*) 2 a + 2 :c*2jta + <■) + (^ - l)«y(c*) 

> (a 2 - 1 ) 4(c*) 2 a + 2c* ^ + a) + (^- - l)c? 2 (c*) 

= (a 2 -l)(2(Q 2 (c*)-c*) + a 3 ) >0 

(note that a,j+ 2 /a.j > a 4 /a 2 ). Therefore, (28) is also satisfied for j + 1, and this completes 
the proof of Theorem 2 . 

5. Some Numerical Examples 

In order to illustrate the range of parameters for which the normalized Chebyshev polyno- 
mials ( 2 ) are not optimal for the approximation problem ( 1 ), we present a few numerical 
examples. Let r* = r*(n) denote the smallest r > 1 such that for all r > r* there exists a 
real number c(r,n) > a r such that for all a r < c < c(r,n) the polynomial ( 2 ) is not best 
possible in ( 1 ). For later use, let us denote by c*(r,n ) the maximal c(r, n) with this prop- 
erty. Recall, that in view of Theorems 1 and 2 , 1 < r*(n) < 00 exists for all integers n > 5. 
In Table I, the numerically computed values of r*(n) and the corresponding semi-axes of 
£ r * are fisted for 5 < n < 20. 


Table I 


Note that r*(n) tends to 1 as n increases. 

The case that the normalized Chebyshev polynomials (2) are not optimal for (1) occurs 
only for c close to the ellipse. In Figure 1, for the cases n — 5 (solid line), n = 7 (dashed), 
n = 10 (dashdot), and n = 15 (dotted), the curves 

c*(r,n) - Or 

(X jr 


10 


are plotted as functions of a r . 


Figure 1 


For some cases for which (2) are not optimal for (1), we computed the best polynomials 
numerically. We were not able to detect any analytic representation of these polynomials. 
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n 

r* 

cz r * 

K* 

n 

r* 

CL r * 

6 r * 

5 

2.6492 

1.5133 

1.1359 

13 

1.3402 

1.0432 

0.2970 

6 

2.0588 

1.2723 

0.7865 

14 

1.3111 

1.0369 

0.2742 

7 

1.8006 

1.1780 

0.6226 

15 

1.2867 

1.0319 

0.2547 

8 

1.6490 

1.1277 

0.5213 

16 

1.2658 

1.0279 

0.2379 

9 

1.5476 

1.0969 

0.4508 

17 

1.2478 

1.0246 

0.2232 

10 

1.4745 

1.0764 

0.3982 

18 

1.2321 

1.0219 

0.2103 

11 

1.4191 

1.0619 

0.3574 

19 

1.2183 

1.0196 

0.1988 

12 

1.3755 

1.0512 

0.3242 

20 

1.2061 

1.0176 

0.1885 


Table I 






